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1  The  purpose  of  «►*«■  work  haa '  been  to  study  the  collective  elementary 
excitations,  as  woll  as  the  ground  state  properties,  of  the  non- re  la t i v is t ic 
many-bedy  system  of  bosons  interacting  with  short  range  forces.  The  ultimate 

/  i  !  >>  1  J. 

aim  is  the  theory  of  superfluid  helium,  a  nd  ,  we—4ws*e  t  he  re  f  o  ro^basM*  concerned 
with  avoiding  the  limitations  to  low  density  and/or  weak  forces  which  charac¬ 
terize  the  Bogoliwbov  theory  of  the  non-ideal  bose  gas  and  those  treatments 
which  derive  from  it.  In  contrast  to  the  fermion  case,  whe-e  the  statistics 
tend  to  suppress  the  many-partic le  correlation  effects  of  short  range  forces, 
there  is  in  the  boson  case  no  reason  to  suppose  that  a  free  particle  or  pair 
approximation  provides  an  adequate  starting  point  for  a  perturbation  procedure, 
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except  in  the  low  density  limit.  dw  basic  approach,  therefore,  has  fcsan  to 
explore  the  use  of  field  theoretic  representations  which  include  from  the 
star^  at  least  some  collective  many-body  effects  of  the  strong  ekor t-range 
interaction. 
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Quaai- Hydrodynamic  R^presentat ion 

The  many-body  Hamiltonian  for  interacting  bosons  (in  the  second- 
quantization  formalism), 
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(  V  i  )*  (Vt)dT  + 
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V(  7  -  r’  )  Yf(  7  )  'ft(  7'  )  Y  (  r')  Y  (  r  )  dT  dT ' 


(1) 


can  be  put  into  a  form  paralleling  the  Hamiltonian  of  a  classical  irrota- 
♦ional  fluid  by  going  over  to  the  action  and  angle  variable  operators  p  , 
S  ,  given  by 

Y  =  e1S/V/i'  ;  Yf  -  p1/2e-1S/a  . 


(2) 
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In  the  c-number  limit  the  transformation  (2)  is  canonical;  however,  as  an 
operator  transformation  it  is  only  spprox  ima  tel  y  unitary-  in  treating  p(  ?  )  , 
St  t*  )  ,  or  their  Fourier  tranforms  , 

cal  commutation  relations 

[  p  (  r  )  ,  S(  r '  )  ]  =  ®  [(  r  -  ?) 

[  p-C,  1  *  I  ‘S,  -k'  (3) 

one  is,  in  fact,  ignoring  the  posi t ive -semidef ini te  character  of  p  .  Thus 
the  representation  is  valid  only  for  states  ir.  which  the  density  fluctuations 
are  small  compared  to  the  mean  density,  i.e. 

<  "k  4  > « o0  (4) 

In  the  case  of  weak  interac tion,  for  states  exhibiting  Bose  conden¬ 
sation,  the  quasi-hydrodynanic  representation  can  be  shown  to  yield  an 
approximation  equivalent  to  the  Bogoliubov  treatment1  o *  the  weakly  non-ideal 
Bose  gas,  the  Bogoliubov  quasi-particles  corresponding  to  density  fluctuations. 

For  the  case  of  strong  Interac  tion  including  a  repulsive  core  of 
radius  a  ,  the  condition  (4)  is  clearly  violated  for  k £  2*/a  ,  since  the 
repulsion  imposes  strong  anti-correlations  in  the  density  over  distances  of 
the  order  of  a  .  However,  at  high  densities  one  can  obtain  an  approximate 

representation  by  working  with  smoothed  density  operators  p  ,  averaged  over 
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volumes  of  the  order  of  r  (r  mg  a)  ,  and  expasiding  the  Hamiltonian  in  powers 
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of  the  Pk  .  .  This  amounts  to  an  expansion  in  Inverse  powers  of  (pQ  rQ) 
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in  contrast  to  the  low  density  expansion  of  Lee -Huang -Yang  for  the  hard-sphere 
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Bose  gas,  which  is  in  positive  powers  of  (p  a  )  ‘  .  In  lowest  (quadratic) 

order  in  the  p,  ,  S,  one  obtains  a  oum  of  harmonic  oscillator  terms  describ- 
k  k 

lng  uncoupled  density  fluctuations  (phonons).  The  excitation  spectrum  is  of 


as  operators  satisfying  canonl- 
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<;no  go r ora  1  form  obtained  by  Feynman'  and  by  Horo.’  Jubov  and  Zubarev, ^  but  diffors 
n  dot  iJ  J  boe-di.-so  of  tho  uwoothing  procedure.  It  ro produced  aumi-quunt 1  tntlvo  J  y 
the  experimentally  obsorvod  excitation  spectrum  for  liquid  Ho  II,  including  tho 

M  t<  8  ‘--l 

roton  minimum  at  tho  obsorvod  wavonumbor  k  k  1,95  x  10  cm 

U 

Tlilu  work  wen  reported  at  tlie  Vllth  Interna tionul  Conforonco  on  l,ow- 
Temperature  PhyslcB,  and  tho  basic  ideas  of  tho  quag i-hyd rod ynamic  represents- 

g 

tlon  are  bIbo  discuased  elsewhere.  --  — 

Even  though  In  lowest  order  the  quasi-hydrodynamlc  representation  gives 
an  adequate  semi-quantitative  description  of  the  collective  elementary  excita¬ 
tions  in  a  Bose  fluid  at  high  densities,  it  does  not  turn  out  to  provide  a  uaable 
starting  point  for  a  perturbation  treatment  applicable  to  the  case  of  liquid 
helium.  The  reason  lies  in  the  largeness  of  the  density  fluctuations  near  tho 
roton  minimum  kQ  which,  in  turn,  is  an  effect  of  the  hard  core  repulsion. 

Only  for  wavenumbers  k  <<:  k^  ,  that  ia  for  distances  much  larger  than  a  ,  is 
a  description  in  terms  of  uncoupled  density  oscillators  an  adequate  starting 
point  for  a  perturbation  treatment.  It  thus  appears  that  an  adequate  treatment 
of  liquid  helium  must  include  essential  density  correlation  effecta  of  the  hard 
core  repulsion  from  the  start. 


Pseudospin  Model 

Considerations  of  the  kind  outlined  above  led  ua  to  investigate  the 
mar.y-body  ground  state  and  elementary  excitations  of  a  simple  "pseudoapin"  cell 
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model,  which  had  been  used  previously  by  Matsubara  and  Matauda  for  treating 
lambda  transition  in  liquid  helium.  Even  though  the  model  haa  some  unphysical 
features — primarily  an  artificial  "band"  structure  due  to  the  use  of  lattice 
quantization,  resulting  in  anisotropy  of  and  a  quasi-momentum  cutoff  in  the 


excitation  spectrum— it  haa  the  advantage  of  allowing  us  to  treat  exactly,  in 
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zero  order,  Important  muny-body  effects  of  the  hard  cor  repulsion,  thus 
overcoming  some  of  the  limitations  of  tho  earlier  treatments. 

The  essential  features  of  tho  pseudospin  model  and  of  our  treatment 
are  those  : 

(1)  Tho  lard  core  ropulsion  is  treated  kinematically,  that  is  b> 

algebraic  relations  txstwoer.  the  field  oporutors,  instead  of  us  the  limit  of  a 

repulsive  potential  in  the  Hamiltonian.  Tho  algebra  of  tho  modol  represents 
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an  approximation  to  Siegert'a  exact  but  highly  untractable  algebra  for  hard 
core  boson  fields.  In  the  cell  approximation  this  algebra  turns  out  to  be 
equivalent  to  that  of  a  set  o  Pauli  spins,  and  the  dynamics  of  tho  model  is 
equivalent  to  that  of  a  Heisenberg  ferromagnet  with  an  interaction  which  is 
anisotropic  in  pseudospin  space.  In  this  model  the  transf ormation  from  field 
amplitudea  to  denaitiea,  correaponding  to  the  action-angle  variable  transfor¬ 
mation  (2),  is  indeed  unitary:  it  ia  aimply  a  rotation  in  psoudospin  space. 
Different  densitiea  pQ  correspond  to  different  orientations  (in  pseudospin 
apace)  of  the  total  "magnetization". 

(2)  (Xir  treatment  splits  the  Hamiltonian  into  two  parts:  an  isotropic 

"unperturbed"  part  Hq  ,  describing  a  system  with  hard  core  ropulsion  plus  a 

2  2 

square-well  attraction  of  depth  h  /md  (d  =  size  of  the  cubical  cell)  between 
nearest' neighbor  cell  pseudOBpins,  and  a  perturbation  consisting  of  the  attrac¬ 
tion's  deviation  from  this  value.  Th®  unperturbed  interaction  corresponds  to 

zero  two-particle  scattering  length  f  ,  and  the  basic  idea  of  our  perturbation 

o 
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treatment  is  thus  similar  to  that  employed  by  Moszkowski  and  Scott  for  the 
nucleon  system. 

(3)  The  pseudospln  algebra  of  the  field  operators  permits  construction, 
by  rotation  of  the  total  pseudospin  from  the  vacuum,  of  r  "unperturbed  q jasi- 
particie  vacuum"  (a  many-body  ground  state  of  Hq)  and  corresponding  quasi-particle 
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oporutors  for  any  moan  density  of  the  Byutom.  Thoso  lorvo  us  starting  point 
for  t routing  tho  full  Hamiltonian  by  tho  oquations-of-motion  mothod  in  the 
random  phaso  approximation.  In  this  approximation  tho  ground  state  of  tho 
full  Hamiltonian  is  obtained  from  tho  quasi-particle  vacuum  (not  tho  fully 
Doso-condonsod  f roo-particlo  ground  state)  by  a  transformation  of  the  Bogoliubov 
typo.  Wo  obtain  expressions  for  tho  excitation  spectrum  and  the  ground  state 
energy : 


Eo/N  -  (2X  ft2  pfQ/m)  ^1  +  (128/15*1/2)(pf3)1/2  (1  -  J^) 
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+  (4n  C  /3)  (pfV/3  ?2/3  (1  -  5  /2) 


} 


(5) 


where  ?o  -  pd  is  the  fractional  depletion  of  the  Bose  condensate  in  tho 

unperturbed  ground  state  (an  "excluded  volusie"  effect  of  the  hard  core 

repulsion),  and  C 2  ■  0,516  is  a  numerical  constant.  In  the  low  density 

limits,  whore  ?o  0  ,  our  results  agree  exactly  with  the  well  known  hard 

2 

sphere  results  of  Lee,  Huang,  and  Yang,  but  the  higher  density  corrections 
reflect  the  raany-body  correlations  due  to  the  hard  core  whioh  are  included  in 
our  unperturbed  ground  state. 


This  work  has  alBo  been  reported. 
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Two-Fluid  Hydrodynamics 

We  have  investigated  the  nature  of  the  velocity  constraint  Introduced 

by  C.  C.  Lin11  into  the  Eckart  variational-principle  derivation  of  two-fluid 

12 

hydrodynamics  for  liquid  He  II,  and  its  effects  on  the  two-fluid  equations 
of  motion  and  the  condition 


curl  v  *  0  . 
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It  appears  that,  while  a  Lin  type  velocity  constraint  must  be  Introduced  if 
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the  Eckart  principle  is  to  be  used  in  connection  with  ordinary  (one-fluid) 

hydrodynamics,  the  form  this  constraint  should  tnko  in  the  two-flu, i  modol 

is  somewhat  arbitrary  and  depends  on  the  basic  assumptions  underlying  the 

model.  In  Lin's  treatment  the  center  oi  mass  velocity  out  not  the  relative 

velocity  of  the  two  fluids  io  constrained.  As  a  result  the  equations  of 

motion  differ  from  those  obtained  earlier  (without  the  Lin  constraint)  and  the 

superfluid  motion  is  no  longer  necessarily  i rrotat ional .  On  the  other  hand 

the  basic  ideas  underlying  the  two-fluid  modol  suggest  that  a  Lin  type  constraint 

should  be  imposed  on  the  normal  fluid  velocity  v^  ,  but  not  on  v^  .  In  this 

case  the  equations  of  motion  are  identical  with  those  obtained  earlier,  and 

curl  v^  remains  zero.  This  work  has  been  reported  in  an  Interim  technical 
13 

report. 

As  a  by-product  of  this  Investigation  Mr  Whitlock  developed  a  varia¬ 
tional  principle  formulation  of  the  lawa  of  motion  for  a  classical  perfect 

14 

compressible  plasma,  which  haa  also  been  reported. 
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